CMB spectra and bispectra calculations: making the flat-sky approximation rigorous 
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This article constructs flat-sky approximations in a controlled way in the context of the cosmic 
microwave background observations for the computation of both spectra and bispectra. For angular 
spectra, it is explicitly shown that there exists a whole family of flat-sky approximations of similar 
accuracy for which the expression and amplitude of next to leading order terms can be explicitly 
computed. It is noted that in this context two limiting cases can be encountered for which the 
expressions can be further simplified. They correspond to cases where either the sources are localized 
in a narrow region (thin-shell approximation) or are slowly varying over a large distance (which leads 
to the so-called Limber approximation). 

Applying this to the calculation of the spectra it is shown that, as long as the late integrated 
Sachs- Wolfe contribution is neglected, the flat-sky approximation at leading order is accurate at 1% 
level for any multipole. 

Generalization of this construction scheme to the bispectra led to the introduction of an alternative 
description of the bispectra for which the flat-sky approximation is well controlled. This is not the 
case for the usual description of the bispectrum in terms of reduced bispectrum for which a flat-sky 
approximation is proposed but the next-to-leading order terms of which remain obscure. 

PACS numbers: 98.80.-k 



I. INTRODUCTION 

Cosmological surveys in general and the cosmic mi- 
crowave background (CMB) in particular are naturally 
constructed on our celestial sphere. Because of the sta- 
tistical isotropy of such observations, cosmological statis- 
tical properties, such as the angular power spectra or the 
bispectra, are better captured in reciprocal space, that is 
in harmonic space. In general however, most of the phys- 
ical mechanisms at play take place at small scale and are 
therefore not expected to affect the whole sky properties. 
For instance, the physics of the CMB is, to a large extent, 
determined by sub-Hubble interactions corrsponding to 
sub-degree scale on our observed sky. Decomposition in 
spherical harmonics, while introducing a lot of compli- 
cation of the calculations, does not carry much physical 
insight into these mechanisms so rather blurs the physics 
at play. 

In this respect, a flat-sky approximation, in which the 
sky is approximated by a 2-dimensional plane tangential 
to the celestial sphere, hence allowing the use of simple 
Cartesian Fourier transforms, drastically simplifies CMB 
computations. Such an approximation is intuitively ex- 



*Electronic address: francis.bcrnardcau@cca.fr 
^Electronic address: cyril.pitrou@port.ac.uk 
t Electronic address: uzan@iap.fr 



pected to be accurate at small scales. So far this approx- 
imation is mostly based on an heuristic correspondence 
between the two sets of harmonic basis (spherical and 
Euclidean) which can be summarized for a scalar valued 
observable by [l|, Q 

e(n) = 5>?mir(ft)-> J d 2 ie(i)e iie . (i) 

In the context of CMB computation, the relations be- 
tween the flat-sky and the full-sky expansions have been 
obtained at leading order in Ref. @. However, its valid- 
ity for the angular power spectrum and the bispectrum 
is not yet understood in the general case and the expres- 
sion and order of magnitude of the next to leading order 
terms are still to be computed. 

The goal of this article is to provide such a system- 
atic construction. In particular, we will show that there 
exists a two-parameter family of flat-sky approximations 
for which well-controlled expansions can be built. That 
allows us to discuss in details their accuracy by perform- 
ing the computation up to next-to-leading order. In the 
new route we propose here, we derive the flat-sky expan- 
sion directly on the 2-point angular correlation function, 
instead of relying on the correspondence ([I} . One of the 
technical key step is to expand the eigenfunctions of the 
spherical Laplacian onto the eigenfunctions of the cylin- 
drical Laplacian in order to relate the (true) spherical 
coordinates on the sky expansion to the cylindrical coor- 
dinates of the flat-sky expansion. This approach proves 
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very powerful since it enables to obtain the full series of 
corrective terms to the flat-sky expansion. 

Once the method has been developed, it can be gen- 
eralized to the polarisation and also to the computation 
of the bispectrum. In this latter case, depending on the 
way one chooses to describe the bispectrum, the exact 
form of the corrective terms has not been obtained but 
we can still provide an approximation whose validity can 
be checked numerically. 

Before we enter the details of our investigations, and 
as the literature can be very confusing regarding the flat- 
sky approximations, let us present the different levels of 
approximations we arc going to use. The reason there 
exist at all a flat-sky approximation is that the phys- 
ical processes at play have a finite angular range. In 
case of the CMB, most of physical processes take place 
at sub-horizon scales and within the last scattering sur- 
face (LSS) (to the exception of the late integrated Sachs- 
Wolfe effect) and therefore within 1 degree scale on the 
sky. Let us denote £o the scale, in harmonic space asso- 
ciated with this angular scale. While using the fiat sky 
approximation, the physical processes will be computed 
in a slightly deformed geometry (changing a conical re- 
gion into a cylindrical one) introducing a priori an error of 
the order of l/£o (actually in 1/£q depending on the type 
of source terms as it will be discussed in details below) . 
Another part of the approximation is related to the pro- 
jection effects which determine the link between physical 
quantities and observables. It introduces another layer 
of approximation of purely geometrical origin. For that 
part the errors behave a priori as l/£ if £ is the scale of 
observation in harmonic space. 

The resulting integrals do not lead to factorizablc prop- 
erties as it is the case for exact computations, while a fac- 
torization property can be recovered taking advantage of 
two possible limiting situations. First, for most of the 
small scale physical processes, one can use the fact that 
the radial extension of the source is much smaller than 
its distance from the observer. It is then possible to per- 
form a thin-shell approximation effectively assuming that 
all sources are at the same distance from the observer. 

Another limit case corresponds to the situation in 
which the source terms are slowly varying and spread 
over a wide range of distances, as e.g. for galaxy dis- 
tribution or wcak-lcnsing field. In this case the sources 
support appears very elongated and it is then possible 
to use the Limber approximation which takes ad- 

vantage of the fact that contributing wave modes in the 
radial direction should be much smaller that the modes 
in the transverse direction (but as such the Limber ap- 
proximation can be used in conjunction of the flat-sky 
approximation or not). 

These different layers of approximations proved to be 
useful to compute efficiently the effects of secondaries 
such as lensing, but also of great help for computing the 
effects of non-linearities at the LSS contributing to the 
bispectrum, either analytically 0, 0] or numerically Q 
(see Ref. [9| to compare to the full-sky expressions) 



as well as for the angular power spectrum; see e.g. 
Ref. [l(| for a review and for the relation between the 
flat-sky and the full-sky expansions in both real and 
harmonic space. We shall thus detail the expressions 
and corrective terms of the flat-sky approximation in 
these two approximations. In particular, we recover the 
result by Ref. @ with a different method in the case of 
the Limber approximation. This is a consistency check 
of our new method. 

First, we consider the computation of the angular 
power spectrum in Section [TT] starting with an example 
of such a construction in order to show explicitly how to 
construct next to leading order terms whose correction is 
found to be of the order of l/£ 2 . We then show that this 
construction is not unique and present the construction of 
a whole family of approximations whose relationship can 
be explicitly uncovered. In Section UTTl we present further 
computation approximations, e.g. the Limber (§ IIII A[) 
and thin-shell (§ IIII Bp approximations. While in Sec- 
tions [IT] and Mil we have assumed, for clarity but also 
because it changes the result only at next-to-lcading or- 
der, that the transfer function was scalar, in Section IIVI 
we provide the general case of the flat-sky approximation 
up to next-to- leading order corrections in l/£ including 
all physical effects. Eventually Section [V] considers the 
case of higher spin quantities to describe the CMB po- 
larization. 

We explore the case of the bispectrum construction 
in Section IVII One issue we encountered here is that 
different equivalent paramctcrizations can be used to de- 
scribe bispectra (amplitudes of bispectra depend on both 
the scale and shape of the triangle formed by three £ 
modes that can be described in different manners). We 
thus present an alternative description of the bispectrum 
for which the flat-sky approximation can be done in a 
controlled way. Although we did not do the calculation 
explicitly, next-to-lcading order terms can be then ob- 
tained. This is not the case for the reduced bispectrum 
for which we could nonetheless propose a general flat-sky 
approximation. Similarly to the case of spectra practical 
computations can then be done in the thin-shell approx- 
imation or the Limber approximation. 

II. POWER SPECTRUM IN THE FLAT-SKY 
LIMIT 

A. General definitions 

In the line of sight approach, the temperature <d(n) 
observed in a direction h is expressed as the sum of all 
emitting sources along the line of sight in direction n, 

9(n) = J dr J ( 27r ) 3 /2 w ( fc > r)<S>(k) cxp(ik.x) , (2) 

where x is the position at distance r and angular posi- 
tion (6, ip) and $ is the primordial gravitational potential 
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from which all initial conditions can be constructed [ll[ . 
w(k, n, r) is a transfer function that depends on both the 
wave-number k and the direction of observation h. In 
particular, it incorporates altogether the visibility func- 
tion, r'e~ r , where r is the optical depth, and the time 
and momentum dependencies of the sources. It can al- 
ways be expanded as 



w(k, h, r) = w jm (k, r)(i) 



3,m 



47T 



2j 



(3) 



where the are the spherical harmonics with azimu- 
tal direction aligned with k. The source multipoles are 
defined using the same conventions as in Refs. [H, [l3| 
except that the multipoles are defined here using the di- 
rection of observation whereas in these references it is 
defined with the direction of propagation (23|. As long 
as we consider only scalar type perturbations in the per- 
turbation theory, the source term will only contain Wj m 
terms with m = 0. For instance, the Doppler term of the 
scalar perturbation introduces a term wio etc. 

In order to focus our attention to the geometrical prop- 
erties of the flat-sky expansion, we first assume for sim- 
plicity that the temperature fluctuations do not depend 
on h and are thus only scalar valued functions. The sta- 
tistical isotropy of the primordial fluctuations then im- 
plies that the fe-dependency reduces to a fc-dependency, 
so that the transfer function is of the form w(k, r) = wqq. 
The general case is postponed to Section IIVI 

The two-point angular correlation function of the tem- 
perature anisotropics, defined by 



m = (e(a)e(fi')>A.A'=cosfl, 

is related to the angular power spectrum Ci by 
C e = 2ir / sm0d0Pi(cos9)£(0), 



(4) 



(5) 



where Pi are the Legendre polynomials of order I. The 3- 
dimensional power spectrum of the primordial potential 
being defined by 



($(fc)$(fc')> = S(k + k')P(k) . 



(6) 



we can easily invert this relation to get a one-parameter 
family of correlation functions £ v (0) as 



dk 



^2 k±dk±w(k, r)w*{k, r')drdr' 



(2tt)2 

x P(k)J [kj_ (r'sinv0 + rsin(l - v)0)} 
x exp [ik r (r cos(l — v)0 — r' cos v0)] (7) 

where Jo is the Bessel function of the first kind of order 
and where the star denotes the complex conjugation. 
The index v refers to the parametrization according to 
the two line-of-sight 



x 
x' 



r[sin((l-v)0)n ±1 cos((l-v)0)] (8) 
r'[s'm(v0)n'j_, cos(v0)] , 



where we have defined the two-dimensional vectors n± = 
(cos tp, sin tp) and n' ± = (cos(p + ip),8m(tp + ip)) with 
ip = 7r so that n' ± = —n±. Note that the relation ([7]) has 
been obtained for a fixed value of tf>, the relative angle 
between n' ± and n±, although it could have been left as 
a free parameter. As we will see in the following, flat- 
sky approximations can indeed be obtained for any fixed 
value of v and ijj provided n and h' are close enough to 
the azimuthal direction. In the last part of III CI we will 
briefly comment on the effect of considering ip. 

Any Fourier mode k can then be decomposed into a 
component k r orthogonal to n± and a component fcj^ 
parallel to n± . Its modulus k is thus to be considered as 
a function of k r and k± since 



To finish, we parameterize k± as 

k = (k± cos /?, k± sin j3, k r ) . 



(9) 



(10) 



Note that the Bessel function in the expression Q arises 
from the integration over the angle between k± and n± , 
i.e. /3 — if. This requires to assume that the transfer 
function depends neither on k = k/k nor on h and is 
thus independent of the angle between k± and n± . Note 
that for simplicity we could have chosen to set tp = 
since only the relative angle between k± and n± matters 
in the derivation of the result. 



B. 



A construction case: the v = case 



Before investigating the full family £ v in the flat-sky 
limit, let us concentrate on the particular case v = 0. 
The flat-sky approximation is obtained as a small angle 
limit, i.e. <C 1 while letting ai = 10 fixed. In order to 
expand the Legendre polynomials P((cos0) in that limit, 
we start from their integral representation as 



Pi (cos ( 



1 



exp [£ log (cos + i sin cos tp)] dip. 



(11) 

In the above mentioned limit, it gives the integral repre- 
sentation of Jo(£9), 

JM6) = - [ exp [i£0 cos tp] dtp . (12) 
Jo 

Furthermore, it allows to obtain the subsequent terms of 
the expansion as 



Pg(cos0) = Jo(ag) 



(13) 



o J iM - ^T J o(a<) + — — 
Z 24 12 a.i 



J\(ott) 



again for a fixed at = 10. The existence of such an ex- 
pansion, and its simplicity, is central in the construction 
we present here. It shows that the eigenf unctions of the 
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Laplacian on the 2-sphere converge toward eigenfunctions 
of the Laplacian on the Euclidean plane. 

This expansion is however not optimal. It was already 
pointed out in Refs. H HH that it can be improved 
by choosing the argument of the Bessel functions to be 
qi = L8 with 



<4 



(14) 



instead of £8. The novel expansion can easily be obtained 
by shifting the argument of the Bessel function in the 
right hand side of the relation fp~3|) . It then reads to 
second order, 



P^(cos( 



Jo(a L ) 

2 r 



■h{a L ) Ji{a L ) 



6a i 



(15) 



The correction term in ~ 8 of the expansion ()13l) has 
indeed disappeared and, as a consequence, the first cor- 
rection to the lowest order of the flat-sky expansion is 
expected to scale as £~ 2 , (in the sense discussed below). 
The accuracy of this mapping is numerically illustrated 
in Fig. [1] for i = 4 and £ = 20, and shown to be bet- 
ter than the percent level. As can be appreciated from 



this expression, from this figure and Fig. [3] later on, this 
change of expansion point is a very important step that 
eventually justifies the use of the flat-sky approximation 
in the context of precision calculations. Finally it is to be 
noticed that in order to use correctly the orthogonality 
relations of the Bessel functions, it is necessary to use the 
variable 



Z = 2 tan (6»/2) 

instead of 8 since Z runs to infinity when 8 — > it 
expansion of Pi{cos8) in function of &l = LZ reads 



P^(cos( 



+ 



2L J 



Jo(a L ) 



(16) 
The 

(17) 



Ji(&l) oil Ji(&l) 
6a l 3 



We are then in position to derive the expression of the 
angular power spectrum in the flat-sky limit. The start- 
ing point is the expression ([7]) of the correlation func- 
tion where the ^-dependent part of the integrand is ex- 
panded as exp [ik r (r — r')] Jo(k±r8) + 0(8 2 ). This can 
then be plugged into Eq. (|5|) , using the expansion (fl"5"|) of 
the Legendre polynomials and integrated over &l, using 
Eqs. (|A7IIA8P to give, up to terms in £~ 2 , 



C 



(u=0) 



— / drdr'dk r k±dk ± c ^ lk ^ r P (k)w(k, r)w*(k, r 1 ) x 



2tt 



(18) 



I ,o(L/r,k ± ) - *l 0>1 (L/r,k±) - ^h,o(L/r,k ± ) + -^h tl (L/r, k±) + ^Ii,i(k±,L/r) - ^I Q ,i(k x ,L/r) 
\2r l bLr 12r 6 4r z 2r 

where the functions I n . p are defined in the appendix [X] Integrating on fcj^, we finally get after integrations by parts 



Cf =°> = i- J drdr'dk. 



cxp[ifc r (r — r')] 



24L 2 



2^± 



P{k)w(k 7 r)w*(k 7 r') , (19) 



r 



where D± = d/d\nk± acts on all the terms on its right, 
and where k 2 = k 2 + L 2 /r 2 . This means that we must 
impose the constraint 



L. 



Higher order corrections will introduce higher order 
derivatives of P{k)w{k,r)w{k 1 r') with respect to 
Since the function on which the derivative acts is only 
a function of k, it is understood in the above expres- 
sion, and in the rest of this paper, that dk ± = (fcj_/fc) dk- 
The evolution of the sources is mainly due to the baryon 
acoustic oscillations and thus dlnw(k,r)/dk ~ ?7lss 
where ?7lss is the mean time of the last scattering sur- 
face. If the initial power spectrum is a power law, 
(9 In P(k)/dk cx 1/k. For super-Hubble scales at the time 
of recombination, the corrective terms are dominated by 



the variations of P(k) whereas on sub-Hubble scales, they 
are dominated by the variation of the transfer function. 
Wc thus have 



(20) d k ln[P(k)w{k,r)w*{k,r')] cx 



1/k if fc7 ?LS s < 1 
77lss if kr) LSS > 1. 



(21) 



The scaling of the corrective term can then be inferred 
easily. If Arxss is the typical width of the last scattering 
surface, the exponential term in the integral on k r im- 
poses the constraint £; r ArLss ^ 1- Now, since the power 
spectrum should be decreasing faster than k~ 2 for large 
k, the dominant contribution to the integral arises from 
k r ^ k±. We thus conclude that k r < min(fcx, ^ r hss) 
and that we always have k ~ k±. In the small scale 
regime, that is for £ > r LSS /?/LSS ~ 50 (r LS s = Vo ~ ? /lss 
being the average comoving distance to the last scatter- 
ing surface) then ki%ss <^ 1- According to Eq. (|2"l"T) . this 
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FIG. 1: Accuracy of the expansions (|13|) in thin lines and (|15l) in thick lines. We plot the difference between Pi(cos9) and its 
approximate forms, at leading order in 8 (i.e. Jo(W), solid lines), first order (long dashed lines), second order (short dashed 
lines) for I = 4 and as a function of 6 (left) and £ = 20 (right). 



means that the first corrective term is, in this regime, of 
order (r/Lss Alss) 2 ~ O(10~ 4 ). In the large scale regime, 
(■ <C ^lss/Aj-lss, and thus k ~ k r rj k±. Since for these 
modes I ^ ''lssA/lss, then fcry LSS <§; 1, and this im- 
plies that the first corrective term scales approximately 
as 1/i 2 . 

Concerning the higher order expansion of Eq. (|19p , we 
see that a corrective term of order 1/L n will involve the 
operator Z?^ / L n where p can be any integer, which gives 
formally a systematic way of organizing the expansion to 
any order. 

To summarize, the precision of the lowest order flat- 
sky approximation in the context of CMB calculations 
is of order l/£ 2 for I < 50 and limited to 10 -4 level for 
£ > 50, i.e. they are of order max(£~ 2 , i^ 2 ) with £ = 50. 
We refer to this first correction as being the correction of 
order l/£ 2 , even if it is limited for small scales because 
of the sub-horizon physics [24| . 

Wc finally note that in some cases the results can ac- 
tually be further improved by choosing the argument of 
the Bessel functions to be ai = L8 = \Jl{l + 1)0, since 
then 



P £ (cos0) = J (a L ) 



( 1L 
2L 



(22) 



M a i) _ 2J i( Q! l) 
3 3a? 



Indeed, this expansion also removes the corrections scal- 
ing as £~ x in Eq. (flT)|) since it also does not contain terms 
linear in 9. Now, the flat-sky constraint reads 



k±r = L= yj £{£+!) 



(23) 



With this choice, it appears that for any transfer function 
w(k,r) that is constant and sharply peaked, i.e. that is 
such as w(k,r) = S(r — rLSs)> and for a scale invariant 
power spectrum, i.e. P(k) = 2ir 2 A 2 k~ 3 , the lowest order 
in the flat-sky expansion takes the form 

^ P ( k ) = 17TT^- ( 24 ) 



This is precisely the result that one would have obtained 
with the exact (or full spherical sky) calculation since 



C, 



dkk 2 P(k) [ji(kr LSS )Y 



2-kA 2 



1) 



27r r £ss 



t(£ + l) 



In other words, with this choice, the leading term of the 
flat-sky expansion is exact. For CMB on large scales, this 
is precisely the case since super- Hubble perturbations are 
frozen, leading to the Sachs- Wolfe plateau. Given that 
the initial power spectrum is expected to be almost scale- 
invariant, we conclude that this expansion of Pi(cos9) is 
the best one in the context of CMB computations and we 
usually adopt it. We can argue that the choice kr = L 
is more compact while the choice kr = L gives a bet- 
ter formula only in the thin-shell approximation. In ap- 
pendix [B] we explain how to change the constraint inside 
the expressions obtained. In general however, one cannot 
state which is the best choice and in the following, unless 
stated otherwise, we will use only L. 



C. Generalized flat-sky expansions 

The derivation of the previous section was limited to 
the case v = but, as it appears clearly from the pa- 
ramctcrizations ([5]), it is actually possible to build a two- 
parameter family of flat-sky approximations, the relation 
between which should be examined. 



1. A one-parameter family of flat- sky approximations 

The aim of this subsection is to explore the conse- 
quence of the use of a more general parameterization 
introducing v as a free parameter. Note that we still keep 
ip = 7T although it could be reintroduced at this stage too. 
Calculations for arbitrary values of ip are however signifi- 
cantly more complicated and do not lead to any improved 
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scheme. Then, again, the particular cases v = or v = 1 
correspond to n z = (0, 0, 1) aligned either with n or n' 
and the case v = 1/2 to n z aligned with (n + n')/\/2. 

We can now follow the same approach as in section lTlBI 
The expansion in powers of Z (instead of 9) can be per- 
formed with the variables A = r — r' and 



R = r v + (1 — v)r. 



(25) 



Plugging the expansion (fTT)) into Eq. (|5|) with the defini- 
tion ([7} of the correlation function, one can perform an 
expansion in Z and A, and then, using the orthonormal- 
ity relations of appendix [X] we can compute the integral 
on eff, in function of the I n ,p(k±, L/R). Eventually the 
result reads 



C. 



(v) 



drdr'dfc r exp(ifc r A) 



2 71 



R 2 



{ ^P(k)w(k,r)w*(k,r'), 
(26) 



where 0^ is an operator that applies on the right part 
of the previous expression with 



O 



24L 2 



[(4 - 3/ 1/2 )i? + /i/ 2 / A] 



D* 



D 



2ARL 2 
D 2 



[/1/2/0A + 3i?/ 2 ]— ^- + ik r (4f R + / 1/2 A)-^(27) 



12RL 2 



where fi/ 2 = 4u(l — v), /o = 1 — 2v and 



8L 2 



D< =fc. 



9 jfc? 9 L 2 d 



dkj 



k dk kR 2 dk 



(28) 



Besides the v = case for which this expression sim- 
plifies, this is also the case for the symmetric choice, 
v = 1/2, for which it leads to, 



C, 



("=1/2) 



drd?''dfc r cxp(ifc r A) 



2n 



R 2 



24L 2 



r 



r n 2 
8L 2 V± 



P{k)w(k,r)w' k (k,r') 



(29) 



The lowest order of this expression matches the one 
derived in Ref. [3]. We remind that the constraint 
k±R = L must be satisfied and that this expression is 
valid only up to order L~ 2 . In the general case for which 
v 7^ 1/2, at lowest order in the flat-sky expansion, the 
expression would remain formally the same, considering 
that R would then be given by R = r'v+(l—v)r. We will 
compare these different parameterizations in the follow- 
ing paragraph. However, if we also take into account the 
corrections, the choice of v would change the expression 
of the flat-sky expansion. In appendix [B] we detail how 
the different corrections obtained are consistent with one 
another. 



In terms of the two-point angular correlation functions 
we have in general, 



a 



dn 1 dn 2 ^(n 1 ,n 2 )F™ 1 (n 1 )y; 2 m2 (n 2 ) ; 



(30) 

where £(ni,n 2 ) is the correlation function of two given 
directions. The difference with Eq. (H} is that we have 
not assumed isotropy. If we choose these two directions 
to be close to a common direction n, then we can expand 
the spherical harmonics of Eq. ()30D around that direction 
which is chosen to be aligned with the azimuthal direc- 
tion and this step breaks explicitly the isotropy of the 
problem. At lowest order, this gives (see appendix |C|) 



2. Breaking of statistical isotropy and off diagonal 
contributions 

The existence of mathematically equivalent flat-sky ap- 
proximations may appear surprising at first view (at least 
it surprised the authors) but it can be fully understood 
when one addresses the construction of the correlators in 
harmonic space not assuming the statistical isotropy of 
the sky. 

Indeed an important consequence of the flat-sky ap- 
proximation is, because a particular direction has been 
singled out, to break the statistical isotropy of the sky. 
There is therefore no reason to have, Ci imi i 27n2 = 



Ce 1 Si 1 e 2 5 ri 



at any order in flat-sky approximation 



where Ct imi e 2 m 2 i s the ensemble average of the product 
of two spherical harmonics coefficients. 



It implies that 



2^+1 

47T 



{-i) m ,j m {tey 



as a consequence of the preservation of the statistical 
rotational invariance around the particular direction h. 
The term Cg^i 2 is explicitely given by 

Ctxi 3 = ^ J dr 1 dr 2 dk^k^dk r P(k)w(k 7 ri)w*(kr 2 ) 

x exp[ifc,.(ri - r 2 )j— j== (32) 

where we remind that L = I + 1/2 (or = + 1)). 
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Performing the integrals on n and T2 leads to 



Ci 



1 

2^ 



dfc 



■Ak r P(k) exp 



i(L, 



^' ¥7 



\fL\L2 

In Fig[2]we present in the space of (£1 + ^)/2 and 

i^l— £2)/%- The main power is carried by multipoles such 
that l\ « £2- The difference between the different possi- 
ble flat-sky approximations lies precisely in the existence 
of (small) off-diagonal terms. 



0.8 



£ 0.4 



0.2 



0.0 



C( t , normalized to unity on ( 




41 K) 



FIG. 2: C ei t 2 /C e+ £ + with 1+ = (Li + L 2 )/2 - 1. Most of 
the signal is localized on the diagonal L\ =1/2- The choice of 
the path of integration in this space leads to different flat-sky 
expansions, that is to different choices of the flat-sky direction 
n or equivalently of the parameter v. For v = or v = 1, the 
lowest order of Eq. (|19[) can be recovered by integrating on 
an horizontal or a vertical line in the {L\,Li) plane. As for 
v = 1/2, the lowest order of the expression (|29|) is recovered 
by integrating on the line L\ + L2 =const. The vertical black 
line(s) is the superposition of those various integration paths 
for L = 280. On the plot they are hardly distinguishable. 



3. Recovering the different flat-sky approximations 

Interestingly, each can be recovered by a proper 
integration of C^ 2 , given by Eq. (f3"2"j) , in the (£1,^2)- 
plane. Each path of integration corresponds to a way to 
relate the correlation function ^(ni,^) to defined 
in Eq. Q. 

For instance, when v = 0, the expression (|19[) . which is 
then the lowest order part of the expression (f?0|) . can be 



recovered by integrating along the path Li=const., i.e. 
as 



LiCe 1 — J &L,2\J LiL 2 Ce 1 i 



(34) 



The integration along L2=const. corresponds to the case 
v = 1 and gives by symmetry the same result, which was 
expected since it corresponds to exchanging h\ and h%. 

The symmetric case v = 1/2 can be recovered by inte- 
grating on the path L% + L2=const. Making the change 
of variables 



La 



we thus obtain 



L\ + L 2 



L_ = 



L\ — L 2 



(35) 




(36) 



Actually, for a general v, the lowest order of the ex- 
pression ([^1) is recovered by integrating on the path of 
equation 



La/Li = (v-l)/v 



(37) 



in this plane. The case of ip ^ tt in Eq. © leads to 
similar constructions but with slightly more complicated 
integration paths. To show it let us reintroduce ip in 
Eq. ((8j) . Then the angular separation 9 of the directions 
n and h' reads 

9 2 = (v 2 + (1 - vf - 2u(l - v) cos ip)9 2 (38) 

and the effective radius R (as it appears in the argument 
of Jo) is 

R 2 = (r 2 v 2 + r' 2 (l -v) 2 - 2rr'v{\ - v) cosvb)6 2 (39) 
with the relation 

L6 = k x R. (40) 

The relation between L, L\ and L2 then follows from the 
relations L± = rk± and L 2 = r'kj_. It reads, 



_ v 2 L\ + (1 - v) 2 L\ - 2i>(l - v) cos ^1^2 



L = 



v 2 + (1 — v) 2 — 2v(l — v) cos-0 



(41) 



which describes the arc of an ellipse in the (L±, L2) plane. 
All these approximations are a priori of similar precision. 
Note however that when tp is small, 9 should be large in 
order to keep 9 fixed making the various expansions less 
precise. It then corresponds to a very squeezed ellipse in 
the (Li, L2) plane. 
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III. TWO LIMIT CASES 

It is worth remarking that in either Eq. (|19j) or more 
generally in Eq. (|29p . the computation of power spectra 
at leading order involves a genuine 3D numerical integra- 
tion. It is therefore numerically less favorable than exact 
calculations (which involves only two integrals [HI, [l|| 
PH-) It is just more transparent since it does not re- 
quire the computation of spherical Bessel functions, and 
geometrically more transparent since the results are pre- 
sented in a Cartesian form. 

Depending on the physical situation, it is however pos- 
sible to introduce simplifications that will make the com- 
putations faster. 



by 2TrS(r — r') so that the integral on r' can then be 
performed trivially. We finally obtain, at lowest order 
in powers of i~ l 



Cp 



dr 



w(k, r) 



P(k), 



(42) 



with the Limber constraint kr = L. The first corrections 
scale as L~ 2 and have been computed in Ref. We 
can recover this result from Eq. (|29[) . by expanding all 
functions of k around k± as 



2k±_ dk 



k—ki 



A. The Limber approximation 

Let us first consider the situation in which the sources 
stretch in a wide range of distances Ar and vary 
smoothly. As long as £/r x Ar ^> 1, it implies that 
k±Ar 3> 1. As a consequence, the contributions in the 
integral on k r are approximately ranging from — 1/Ar to 
1 / Ar , that is in a range of values where k = \J k 'j_ + fc 2 . ~ 
k±_ and all functions except the exponential in the inte- 
grand can be considered constant. 

The Limber approximation then consists in replacing 

J dk r exp[ifc r (r — r')] 



This will result in a term proportional to k 2 which can 
be handled using 

J dfc r (i£v)"exp[ifc r A] = 27T<S { ' l} (A) , (43) 

where 5^ is the n-th derivative of the Dirac distribution. 
The term proportional to k r in Eq. ([29)) can be handled 
in the same way and gives 8' (A). Integrating by parts 
in A removes the derivatives of the Dirac distributions, 
using that dR/dr = 1 — v and dR/dr' = v. Then, the 
integral on r' is trivial, because of the Dirac distributions, 
and one can then perform an integration by parts in r in 
order to reshape the result. We finally obtain 



C, 



(•"=1/2) 



dr 



D 



-} P (k)\w{k,rW + — 2 



d[w(k,r)/y/r\ 



<91nr 



P(k) 



^ [w(k, r) 2 P{k)] - ^ [w(k, r) 2 P{k)] 1 (44) 



where D n = k n d/dk n and we recover the result of 
Ref. [f|. The Limber approximation is not well suited 
for CMB computation since its hypothesis is not satis- 
fied during recombination. However, when it comes to 
the contributions of the reionization era, sources stretch 
in a wide range of distances and the Limber approxi- 
mation can be used. Special care must be taken though 
because the source have a directional dependence and are 
not simple scalar valued sources. The main sources are 
(i) the Sachs-Wolfe contribution <5 r /4 + $ where 8 r is the 
density contrast of the radiation (it) the late variation 
of the gravitational potential, and (Hi) the Doppler con- 
tribution n z diV, where v is the scalar part of the baryon 
velocity. A naive estimation would lead to think that the 
Doppler effect is the most significant contribution since 
v ~ (rLss — r ) while the intrinsic Sachs- Wolfe term re- 
mains of order one. However, in Fourier space this leads 
to a term proportional to k r v, and it cannot be dealt 
with the naive replacement k r — > since this would van- 



ish. Instead, it should be treated using Eq. (|4"3")l and 
it leads to replace the Doppler source by idv/dr (if we 
ignore the derivatives of 1/r) before taking the lowest 
order of the Limber approximation, and this then gives 
a contribution of order unity as well. Physically, this 
means that the modes which favor the Doppler effect arc 
those aligned with the line-of-sight n 1 , but the contribu- 
tion of these modes is further suppressed in the Limber 
approximation for sources stretching in a wide range of 
distances. See Section IIVI below for a discussion on the 
flat-sky expansion in general in the case of sources with 
intrinsic directional dependence. Furthermore, if we de- 
cide to consider the late Integrated Sachs- Wolfe (ISW) 
effect which is due to the variation of the gravitational 
potentials, the Limber approximation is satisfactory as 
well. We shall not consider any of these in this article 
and focus our attention to the effects that occur during 
the recombination. 
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B. The thin-shell approximation 



with 



A useful approximation can be derived when the 
sources are contributing only in a thin range of distances. 

Indeed, in such a case the factors 1/r 2 or 1/R 2 can be 
replaced by l/?Lsg, wnere r LSS is the mean distance of 
the sources each time one has to compute k (e.g. k = 
^k 2 +L 2 /R 2 is replaced by fc L ss = ^/ k 2 + L 2 / r 2 ss ) . 
The integrals on r and r' in the lowest order term of 
expression can be factorized and we obtain the thin- 
shell flat-sky approximation 



SO 



1 



' LSS 



dk r 
~27 



drexp(ifc r .r)w(fc LSS ,r) 



P(k 



LSS) 



(45) 

In the full-sky computation, this factorization is auto- 
matic since there is also an integral over r, which is then 
squared, and another integral over k. In general, the flat- 
sky expansion breaks this property, and it is recovered in 
the thin-shell approximation. 

Estimating the error introduced by the thin-shell ap- 
proximation is actually difficult and quite model de- 
pendent. First the error introduced by replacing the 
factor 1/R 2 by l/r^gg will be very limited if tlss is 
taken in the middle of the LSS and should be much 
less than ArLss/ r LSS- However we also need to esti- 
mate the error introduced by the approximation of the 
flat-sky constraint. Depending on the k dependence of 
S(k,r) = w(k,r)w(k,r')* P(k), this will lead to two ex- 
treme cases. If it is a pure power law and depends only 
on k, then S(k) ~ S(k hS s) - (R- rhss)/rhssD±S(k hS s) 
and again if tlss is chosen in the middle of the LSS, af- 
ter integration on r and r' , the error would be much 
smaller than ArLss/ r LSS- However if the sources are 
purely oscillatory with frequency k, which is the case 
for the small scales, that is if S(k, r) cx cos(fcr) for in- 
stance, the relative error introduced would be of order 
ArLSs/^LSS- So for small scales, this is much larger than 
the first corrections in the expression (fT9| which scale 



as max(f"Vo ) with £ = Tlss/^lss — 50. We shall 
see in the next section that it is in principle comparable 
on small scales to the larger corrections introduced when 
we consider the non-scalar nature of the sources which 
are of order max(!"Vo ^ However, in practice the er- 
ror introduced in the thin-shell approximation is smaller 
since At^lss < Vlss and also because the sources are not 
purely oscillatory and converge to a power-law on small 
scales thanks to viscous effects, thus reducing further the 
error made. 

We can also try to compute corrective terms to the 
thin-shell approximation. This correction is obtained by 
expanding R around Tlss but also k around /clss and is 
given at lowest order of this expansion by 



^^fthiii shell 

f drdr'dk r exp(ifc r A) 



2tt 



(46) 

SOP(k LS s)w(k LSS ,r)w* (k LSS ,r') , 



(R- r LSS ) 
rLSS 



(47) 



LSS 



Since R-r LSS = v[r' - r LS s] +(l—v)[r — r hss ] the inte- 
grals on r and r' in the expression of this correction can 
also be expressed as a sum of factorized integrals, which 
means that numerically it corresponds effectively to sums 
of two-dimensional integrals. On large scales the operator 
D±_ acts mainly on P(fcLSs) since it contains the domi- 
nant k dependence (sec the discussion in section III Bl) . 
However on small scales the dependence in k is domi- 
nated by the source terms w(fcLsSi r ) and w*(khsST r ')i 
and this lowest order correction is not valid given the nu- 
merous oscillations of the sources. However for a source 
which depends purely on fcr, as is roughly the case for 
the baryon acoustic oscillations, then on small scales it 
depends nearly on k±r since k ~ k±, and the error intro- 
duced by replacing k with /clss can be seen as an error 
in the placement of the distance r at which the source 
is emitting. In that limit case, everything happens as if 
the visibility function contained in the expression of the 
source w{k,r) was slighlty distorted when performing the 
thin-shell approximation. In practice, we shall not cor- 
rect for this since the source is not purely depending on 
kr on small scales. This means that we should make the 
operator D± contained in SO act only on -P(/clss) for our 
practical purposes when using Eq. (|46[) to correct for the 
thin-shell approximation that we take. 

In the context of CMB, the use of the thin-shell approx- 
imation requires to rewrite the source in order to localize 
the physical effects on the LSS. In practice, the terms 
involving the gravitational potential $, i.e. of the type 
n i 9i$ [111, [lH, whose contribution would stretch from 
the LSS up to now, are replaced by d$/dr — <9$/<9r. 
This clearly splits the effect into an effect on the LSS 
(d<I>/dr), the Einstein effect, and an integrated effect 
id^/dr) which is negligible for a matter dominated uni- 
verse since the potential is then constant. Should we 
consider the effect of the cosmological constant on the 
variation of the gravitational potential, then we could 
use the Limber approximation discussed in the previous 
paragraph, and sum the resulting Ci to the contribution 
of the LSS. 

In Fig. [31 we present the flat-sky approximation for 
an instantaneous recombination, including only the in- 
trinsic Sachs- Wolfe effect in order for the source to be 
purely scalar. The corrections in the thin-shell approx- 
imation can be read from the expression (|29[) . However 
for practical purposes, the derivatives with respect to fcj^ 
need to be converted into derivatives with respect to k so 
that their action on w(k, r)w(k, r')*P(k) is clearer. How- 
ever, as we shall see in Section HVl for realistic purposes, 
the sources are not a pure scalar and have an intrinsic 
geometric dependence so that the first correction will ac- 
tually scale as L^ 1 . 

It is interesting to remark that all flat-sky expansions 
will lead to the same expressions at lowest order in the 
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FIG. 3: Left: Comparison of the flat-sky approximations (with k±r = £ in dashed line and with k±r = \J £{l + 1) in dotted line) 
to the exact computation (solid line) . We consider the standard cosmology with an instantaneous recombination and ignore all 
effects but the intrinsic Sachs- Wolfe effect (Q = <5 r /4 + $). Right: The relative errors of these two flat-sky approximations with 
respect to the exact computation. The first method is limited to a 1% relative error above I = 100, as discussed in Section Til Bl 
whereas the second one is much better since the first corrections scale as £~ 2 . 



thin-shell approximation and in the Limber approxima- 
tion, since then, as we have already seen, the lowest order 
does not depend on the choice of v (nor i/j). In the case 
of the thin-shell this can be understood from Fig. [2] The 
thincr the shell, the more the function Gi x ^ is peaked on 
the diagonal L\ = and the less the integral depends 
on the line of integration in the (L\, L2) plane. On Fig. [2j 
we plot two contours of integration corresponding respec- 
tively to v = and v = 1, and it can be understood that 
the integrals obtained on these contours cannot be very 
different. The Limber case is similar since the factor 



dk r exp 



in Eq. 



is approximated to be 27rfcj_i5(Li — L2). 



on small scales, where we remind that f3 is defined in the 
parameterization (|10[) of k. 



A. Lowest order flat-sky expansion 



As long as we consider only scalar perturbations, the 
source term will contain Wj m terms with m = 0. Note 
that this is different from our previous assumption that 
only woo was not vanishing. 

The expansion ([3]) contains only terms in Y£° (n) which 
(48) are proportional to Pj(k.n). At lowest order in 0, 
Pj n (k.h) depends only on k r /k so that the flat-sky ex- 
pansion remains unchanged as long as we make the re- 
placement 



IV. EFFECT OF NON-SCALAR SOURCE 
TERMS 



w(k,r) —> i J Wjo(fc, r)Pj (k r /k) 



(50) 



So far, we have assumed that the transer function was 
scalar, in the sense that it was a function w(k,r) that 
does not depend on k and n, that is the expansion J3]) 
contained only a term u>oo- In general, this is not the 
case since scalar perturbations involve m = terms with 
£ = 1 (Dopplcr effect) while vector perturbations and 
gravity waves generates terms with m = 1 and ni = 2 
respectively. 

In order to take this dependence into account in the 
flat-sky analysis, one needs to compute Eq. (J7J) with the 
source . From the parameterization §E§ with the choice 
tp = 0, wc deduce that 



k.n = i 

k 



k r + k ± -cos(P)+O(0 2 ) 



(49) 



In conclusion, the formal expression of the fiat-sky ap- 
proximation at lowest order remains unchanged for cos- 
mological scalar perturbations (which are the dominant 
sources of CMB anisotropics). 



B. First correction 

We have seen that for scalar field sources, the first 
correction scale as l~ 2 . However, for the general case, the 
first correction arises from the first correction in (|49[) and 
scales as l~ Y . Since the dominant effects come from j = 
and j = 1, we can ignore the contribution coming from 
j = 2 in the computation and there is no contribution 
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for j > 3. Thus, we restrict to 

k 

w(k,h,r) = [w ao (k,r) + i-?-w 10 (k,r)} 



i—wio(k, r) cos /? 



C(6» 2 ).(51) 



In the two-point correlation function, the first correction 
will arise from the product of the first order of w(k, n, r) 
with its lowest order term. Using the symmetry of the 
integral in k r the first corrective term due to the geometry 
of the sources is thus 



w(k, n, r)w*(k, h' , r') ~ (52) 
k 

w o(k,r) + i-^w 10 (k,r) 
9 k± 

+ -icos/3— [w o(k,r)wl (k,r') + w w {k,r)w^ (k,r')] 
= ww*(k, r, r') + 9\ cos (3 {ww k y > (k, r, r') . 



In order to go from Eq. @ to Eq. ([7]), the integral over 
P will give a term — Ji[k±RZ] instead of the previous 
term Jo[k±RZ] (due to the factor icos/3). Following the 
exactly same method as in Section III CI we obtain the 
flat-sky expression with the first correction included 



C, 



drdr'dfc r exp[i/c r (r — r')\ 



2vr 



D± 
k\r 



^P(k)ww*(k, r, r') 
k ± P(k) (ww*) {1) (k,r,r') }, (53) 



where we can take either the fiat-sky constraint k±R = L 
or k±R = L given that the corrections to this expression 
are of order f 2 . This expression can be further sim- 
plified in the thin-shell approximation and can then be 
further improved by including the corrections due to this 
thin-shell approximation given in Eq. (|46[) . Similarly to 
our discussion at the end of Section III Bt the first cor- 
rective term which comes from the second term in the 
curly brackets in the expression above, is of order 1/ £ up 
to £ ~ 100 and then of the order of 1% beyond. Includ- 
ing this first correction to the lowest order of the flat-sky 
expansion (the first term in the curly brackets in the ex- 
pression above) considerably improves the precision of 
the fiat-sky expansion. Note that the corrective terms 
in the expression (|2T))) were at least of order l/£ 2 , and 
thus the corrective term in the expression (|53|) . which 
comes from the directional dependence of the sources, is 
the dominant one. To compare with Eq. (|29[) . one would 
need to derive the corrections of order l/£ 2 in Eq. (f53")l . 
which seems an unneeded academic sophistication. Fig.0] 
depicts the relative error with respect to the exact calcu- 
lation, with and without including the first correction of 
order l/£ given in Eq. (|53|) . using the thin-shell approx- 
imation and taking into account all effects but the late 
integrated Sachs- Wolfe effect. 



V. FLAT-SKY EXPANSION OF HIGHER SPIN 
QUANTITIES 

The CMB radiation is not described entirely by its 
temperature, since it is polarized by the Compton scat- 
tering of photons on free electrons. Only the linear po- 
larization is generated through this process and it is de- 
scribed by the spin ±2 fields defined from the Stokes 
parameters 



I 2 



X{n) ee Q(n)±iU(n). 



(54) 



±2X{n) is dependent on the choice of the basis used to 
defined the linear polarization. Any rotation of this basis 
by an angle 7 around the direction h transforms it as 
±2 X(n) ^ ±2 X{h)e ±2 ^. 

In order to define the correlation function of two 
spinned quantities, Sl X and S2 Y, one shall use the spin 
raising and lowering operators, respectively d + and d~ 
which are defined for a spin s quantity by U 53] 



X = - sin ±s (d 6 ± i esc 0d v ) sin TS S X 



(55) 



in order to relate them to a spin-0 field. We thus define 

s x = (-iyd- s s x, s x = d- s s x (56) 

respectively for s > and s < 0. The expansion of S X 
on spinned spherical harmonics according to 



s X(n) = ^ s X £ms y/ n (n) 



(57) 



/ in 



can then be related to its expansion on spherical harmon- 
ics as 



,X{n)=Y J sX l 



(n, 



Yt(n) 



(58) 



The two sets of spherical harmonics are related by 

d-' B Y? = 



\£_ 



: 1 1> 



if s > 0, (59) 



if s < 0. 



We further define the electric and magnetic parts as 
E(X)(h) = \ [ s X(h) + _ s X(n) 



B{X)(n) = - [ s X(n) - _ s X(n) 



(60) 
(61) 



where here we choose the convention s > 0. For a spin 
quantity such as the temperature, E{ X) — X = X 
and B(qX) = 0. Due to parity invariance, the correlation 
between an electric type quantity and a magnetic type 
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FIG. 4: Left: Comparison of the fiat-sky approximations with the constraint k±r = \J i{l + 1) while including only the lowest 
order in the expression (|53[) (dashed line) or adding its correction of order f (dotted line) to the exact computation (solid 
line, standard cosmology, ignoring only the late time ISW effect). The calculations in the flat-sky approximations have been 
done using the thin-shell approximation. 

Right: The relative errors of these two approximations with respect to the exact computation in respectively continuous for 
the lowest order expression and dashed line when including the first correction in l~ . In dotted line we also compute the error 
when adding on top of the order I -1 correction, the correction of Eq. (|46[) due to the thin-shell approximation. The lowest 
order is limited to a 1% relative error beyond I = 100 whereas the first correction increases substantially the presicion on small 
scales and the correction for the thin-shell approximation improves also the largest scales. 



quantity always vanishes. We then define the correlation 
function as 

£ E (x)e(y) {9) = ( E (X)(h)E(Y){n')*)fi.w=cose (62) 

21+1 j (£ + Si y. (£ + Izjl C E(X)E(Y) 



4tt V (£-si)! (i-s 2 )\ 



>Pi(cos6), 



with similar definitions for magnetic type multipoles. 
The angular power spectra are then extracted through 



r E{x)E{Y) _ l (£- Sl )\(e- S2 )\ 

y (l+srV-il + szy. 
x / sin6»de» ^(cos60£ B ( A ") B ^)(60 



(63) 



The emitting sources for a spin s quantity are expanded 
similarly to Eq. ([3]) but with a decomposition on spinncd 
spherical harmonics 



w[ s X}(x) 



d 3 k 

(2tt) 3 / 2 



w[ s X](k,ft,r) exp(ik.x) (64) 



with 



w 



-in 



,X\(k,n,r) = 5>, ro [_ s xp,r)(i^ . 

£,m 

x a Y£ m (n)er*t h ' A \ (65) 



where (p(k, ft) is the azimuthal angle of k with respect to 
ft [l[. The source multipoles for spinned quantities are 



defined using the same conventions as in Refs. [12J, [l3j 
except that the multipoles here refer to the direction of 
observation whereas in these references it refers to the di- 
rection of propagation. We have thus made the replace- 
ment s — > — s additionally to the extra factor (— 1)^ which 
was already considered for spin quantities in Eq. 
in order to take this fact into account by using the trans- 
formation properties under parity of the multipoles. 

At the lowest order of the flat-sky expansion, we can 
approximate <p(k, ft) by /3. In particular, this implies 
that the spin raising and lowering operators applied on 
u>[ s X] will act only on exp(ifc.a;) and we find that at 
lowest order in the flat-sky expansion 



d exp(ik.x) ~ (-ifcr) s e 



exp(ife.a;), (66) 



for s > 0. We thus deduce that the sources for E( S X) 
and B( S X) are given by 



{(E/B)( s X)}(k,h,r)=Y,(ikry(i) 



l.m 



4tt 
2£+l 



(67) 



±-w im [(E/B)(X)](k,r) 



x [(-l) s e- isv _ s y^ m (n) + e 1 ^ s Yi m {ft)} 
+ ^w em [(B/E)(X)](k,r) 

x [(-i)' e - ia *_ a y^ m (n) - e is ^Yt l (n)} } 
with the + sign for E and the — sign for B. 
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In the case where there are only scalar sources (m = 0), 
this simplifies substantially. From the parameterization 
(jHJ, that is with the choice <p = and using that s Y eo = 
( — l) s Y ls , we obtain 



w[E(X)](k,n,r) 



(68) 



^(-ifc r )«(i)^4S(X)](/c,r)y| T |}p/(fc r /fc) 

and a vanishing w[£?(X)]. Following the same method 
as in section III Bl the factors (kr)( Sl+S2 ' ) are going to 
be approximately canceled by the prefactor of Eq. (fM)) 
which is behaving as l~^ 1+S 2 ). A mo re careful derivation 
would require to use k = \J fc^ + fc 2 . and the exact form 
of the prefactor in Eq. (|64[) . However, we are here inter- 
ested in the most simple expression for the lowest order 
of the flat-sky expansion and we drop these extra com- 
plications. Since polarization is not generated on large 
scales, this will be sufficient to obtain an excellent flat- 
sky expansion. In the end of the computation, for scalar 
perturbations the multipolcs associated with the corre- 
lation of spinned quantities is obtained just by replacing 
the sources according to 



w(k,r) -^^2(-l)'i t+ 'wto[E\(k,r) 



{£ + s)\ 1 \ k 



We use this expression for the computation of the sources 
of linear polarization (s = 2, and only £ = 2 in the sum 
above since there are only quadrupolar sources), and wc 
compare the full-sky computation of Cf E and Cf E with 
the flat-sky result in Fig. [5] 



VI. THE BISPECTRUM 

We shall now investigate the flat-sky expansion of the 
bispectrum. We shall assume that it arises from pri- 
mordial non-Gaussian initial conditions described by a 
primordial bispectrum for the metric fluctuation, e.g., 

(*(fci)*(fc 2 )*(fc 3 )) - <5 3 (fci + k 2 + fc 3 ) B(ki, fe 2 ) • (69) 

As we shall see the derivation of the temperature bis- 
pectrum is much more subtle and we only partially suc- 
ceeded in a sense we explain below. The reason lies in the 
fact that there are actually many ways of representing a 
bispectrum the properties of which in the flat-sky limit 
might be different. 



Two different representations of the bispectrum 
in harmonic space 

1. Definitions 

We are interested in the expectation values of 
i(imiflf 2 m 2 (1 <3m3) °f the a,£ m coefficients of the CMB 



temperature. Because of the statistical isotropy of the 
sky, their m-dependence is bound to be that of the Gaunt 
integral Qmlmlm^ so that it 18 more fruitful to introduce 
the reduced bispectrum b^i 
defined as 

(^imi Qt 2 rn 2 ^37713 I 



(see Ref. [18j for instance) 



^mim2m3 



(70) 



The a priori purpose of the following is to propose a con- 
trolled approximation for bi^i^ in the flat-sky limit. It 
turns out however, as we shall see later, that in order 
to have a controlled limit expression of that quantity, 
strong regularity conditions should be imposed on the 
initial metric perturbation B(ki,k2). 

The bispectrum can actually be equally characterized 
by the following quantities (reminding L = £ + 1/2), 



h £2 £ 3 \(£i £2 t 3 
3 1 M -M H I ' 



(71) 



that contains the same information as fr^&g^ since it can 



be inverted as 



h £2 £3 
K M —M 7 
2 , ~~Tn — 7, — rr&feW' ( 72 ) 



M 



£1 £2 £3 

000 



&i^ 2 m is a real- valued quantity that obeys the following 
symmetry properties 



(73) 



As we shall see, £^ 2 m actually enjoys a better-controlled 
asymptotic expression than bi^l^ m the large £ limit, 
similar to what was achieved for the angular power spec- 
trum in the previous sections. 



2. Properties of fy^M 

Let us first relate Cti£ 2 M to the angular three-point 
function. We introduce three unit vectors on the celestial 
sphere, hi, h% and 77,3, by 



= (sin 9i cos ifi, sin#i simpi, cos 61), (74) 

where 9i and ipi are the Eulcr angles and i — 1..3. Tak- 
ing advantage of the statistical isotropy of the sky, the 
three-point temperature correlation function can always 
be expressed as a function of the relative angle with re- 
spect to say tt-3, setting #3 = 0, so as a function of the 
four angles, (6±, (fi, 62 , f-z)- These dependencies can then 
be expanded on spherical harmonics as 



(6(n 1 )e(n 2 )e(n 3 )) 

= y ] £,e 1 £ 2 m 1 m 2 
(,,m, 



= ^1,^1^2,^2) (75) 

or or 
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t ( 

FIG. 5: Comparison of the fiat-sky approximation (dashed line) to the exact computation (solid line) for Cf E (left) and Cf E 
(right) ignoring the effect of reionization. Since they only disagree on very large scales, in the regime where the polarization 
fails to be generated, both curves are hardly distinguishible. We assume standard cosmology. 



Rotational invariancc further implies that £ depends only 
on (/?2i = ¥>2 — so that £i' £> m ' m i vanishes if m[ ^ —m' 2 . 
We conclude that 



£(0i, Vl, 02,^2) 



(76) 



This expression generalizes the expansion of the two- 
point correlation function in terms of Pe{0). 



3. Relation between ( > i 1 e 2 M and the bispectrum 

To obtain such a relation, we need to express £ in terms 
of directions hi, h 2 , h% instead of relative angles. This 
can be achieved by performing a rotation 1Z, under which 
the sperical harmonics transform as 



Yr(n- 1 h) = Y,Y e m '(h)Di, m (n) 



(77) 



where D m , m are the rotation matrices and can be ex- 
pressed in terms of spin- weighted spherical harmonics as 

/ 47T 

D e _ ms (<p, 9, </>) = (-1)'" y ^ sY e m (d, <p) e"^. (78) 
It follows that 

£(ni,n 2 ,n 3 ) = E &iI*m-m (79) 

^1 ,&2 )M ,1711 ,7712 

xyr*(«i) MY e ^(h 3 )Y^*(h 2 ) - M Y^(h 3 ). 



Using 



d 2 h Sl Y^{h) S2 Y^(h) 83 Y e ^(h) = (80) 



8L!L 2 L 3 f h £ 2 £ 3 \ f h £ 2 I 



Ait \s\ s 2 S3 J \mi m 2 m 3 
when si + s 2 + S3 = 0, we obtain 



h £ 2 £ 3 

V M — M 
2 4 — — — 7, — ~ — £e 1 e 2 M-M- (81) 



h £2 £3 




This shows that the coefficients ^t^M-M ar< 3 nothing 
but the parameters £,t^i 2 M that we introduced earlier as 
an alternative description of the bispectrum, 



iiiliM — £tie 2 M-AI ■ 



(82) 



They can therefore be expressed in terms of the real space 
correlation function, 



47T 



d 2 h 1 d 2 h 2 £(0 1 ,ipi,0 2 ,tp 2 ) 

(83) 

sin#id#i / sin#2d02 / d</? 2 i 



V2Li2L 2 

xYg (e l ,<p 1 )Y l - M (e 2 ,< P2 ) 

Sir 2 



V2ii2L 2 

xYg (0i,o) r- M (^ 2 , ^ 21 ) e(0i,o, e 2 , <p 2l ). 

This is the generalisation of Eq. ([5]) for the bispectrum. 

B. Flat-sky limit of S,i 1 i 2 m 

We follow the same path as for the power spectrum. 
The first step is then to provide a formal expression of 
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the three-point correlation function in real space. Setting 

Xi = nfk, (84) 

for i = 1...3 and 

k\ = (fei sin ai cos/?i, k± sinai sin/3i, k% cosai) 

= (^cosjSi,/ci-sin^i,A;f), (85) 
k 2 = (fc^cos/3 2 ,fc^sin/3 2 ,fc2) ■ ( 86 ) 



Formally the three-point temperature correlation func- 
tion then reads 



£_(9 1} ip 1 ,9 2 ,ip 2 ) = J d 3 fcid 3 k 2 dri dr 2 dr 3 B(ki, k 2 ) exp [ikf(ri cos6q - r 3 ) + ifcf(r 2 cos# 2 - r 3 )] 

xw(ki, n) u>(fc 2 , r 2 )w(|fc 3 |, r 3 ) exp n sin#i cos(/3i — <pi) + ifc^^ sin# 2 cos(/3 2 — <^ 2 )] , (87) 

I 



where the Dirac distribution of Eq. ([69]) has been taken 
into account. We are left with a function of the relative 
angles (9 1 ,<p 1 ) and {9 2 ,p 2 ). 

For a fixed value of M, Y^' I {9, tp) has a well con- 
trolled limit in the flat-sky approximation. It is given 
by Eq. (8.722) of Ref. 



Y, 



Mi 



1/2 



(-l) M JM[£% iMip . 



In appendix [Cj we show how the next to leading order 
terms of this expression can be obtained. The expan- 
sion parameter is M/£ or M9. The completion of the 
calculation then relies on the relation, 



xJ M {ax)J M (bx)dx = ^ . (89) 



We can now proceed to the evaluation of (|83|) in the 
flat-sky limit, as it is now straightforward. Defin- 
ing pi = k^~ri9i and p 2 = fc 2 r 2 ^ 2 the expression of 
<pi, 62, ^2) at leading order is 



£(01,0,02,^21) = 

J d 3 fc!d 3 k 2 dn dr 2 dr 3 B{k l , k 2 



(90) 



exp 



i V k 



2^ K a r a 3 
0=1,2 



w(k a ,r a 

=1,2,3 

exp [ipi cos/?i + \p 2 cos(/3 2 — <y9 2 i) 



with rj 3 = r, — r 3 , and fc 3 ee |fe 3 |. Then the subse- 
quent angular integrations that appear in the expression 
of &i£ 2 m in Eq. (|83|) lead to the following transforms: 

• the integration over tp 2 i of 
exp [ip 2 cos(/3 2 — — \Mp 2 {\ gives a term 
i M J M (p 2 )e iM *; 

• the integration over (3i, fixing the relative 
angle j3\ 2 between the wave vectors k\ and 



k 2 , of exp [ipi cos fii + iMf3 2 ] gives a term 
(-i) M JM(pi)e iM fe; 

• the integration over 6q of Jm{pi)Jm{L]_9i) gives a 
term 8{L\ — k^rij/Lx] 

• the integration over 2 of Jm{P2)Jm{L 2 9 2 ) gives a 
term <5(L 2 — k 2 r 2 )/L 2 . 

Then, the integration over k^ and k^ can be performed 
explicitly and we are left with 



where 



„*, , , w(fci,ri) w(k 2 ,r 2 ) 

dk 2 dri dr 2 dr 3 2 2 

r i r 2 
xw(fc 3 , r 3 )B(k u k 2 ) exp [ifc*ri 3 + ifc 2 r 23 ] (92) 



U 3 (/>{*}) = / d*J 



In this expression, 



fcf + L*/r£ (93) 



and B(ki, k 2 ) is an implicit relation of fef, fcf and of the 
relative angle of their transverse parts, since 



kl = {kt + klf + Ll/rl+Ll/rl 
+2cos(/3 {k} )L 1 L 2 /(r 1 r 2 ). 



(94) 



It can further be noted that because &f^ 2 (/3{fc}) is invari- 
ant under /3rj.i — > — the expression of £^ 2 m also 
reads 



d/3 {fe} 



2^ 



3 (M/3 {fc} ) (%- } ). (95) 



The M-dependence in this expression is the one of the 
Fourier transform of the /3{^}-dependence, that is that 
of the relative angle between the wave numbers in the 
transverse direction. The equations (l§2"j) and (jHSj) repre- 
sent the flat-sky approximation of ^t^M- 
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C. The flat-sky limit of be 1 e 2 e 3 

While it was straightforward to derive the flat-sky limit 
of ^,< 2 m, the one of b^ 2 ^ 3 is more problematic. First, 
biii 2 i 3 is obtained from a sum of contributions each of 
which involves ^i 1 i 2 m which should be calculated in the 
flat-sky limit. If the number of M in that sum is finite, it 
is however possible to obtain bi 1 i 2 i a . There is however a 
priori no reasons for the sum in Eq. (|8ip to be dominated 
by its first terms. It depends actually on the details of the 
model and in particular on the regularity of B{ki 1 k,2). 
Fortunately, this should be the case for the models of 
interest in cosmology. 

We can then try to invert Eq. (|7Tj) in the large £ 3 limit . 
First, Eqs. (6.578.8) and (8.754.2) of Ref. [H (which are 
a particular case of the Ponzano and Regge semiclassical 
limit of the Wigner coefficients [2(|) allow to infer the 
limit 



This means that it fixes = /3{l} an d thus fixes the 
value of k^ as a function of r 3 and L 3 . This relation is 
actually nontrivial as it is not given by k 3 = L 3 /r 3 as 
one would naively expect. It is instead given by 



fi, 3 — 



Li (ri 



rir-2 r\T2 \r\ 



— - 1 



7-2 

L 2 



n-i 



(101) 



If the last scattering is thin enough so that we can ap- 
proximate r\ = T2 = r 3 , that is in the thin shell approxi- 
mation, we recover what we would have naively expected, 
i.e. k 3 L r 3 = L 3 . 

Similarly to the angular power spectrum case, we have 
used L = £ + 1 /2 instead of £ to improve the convergence 
of the flat-sky expansion. We have not been able however 
to infer formally the validity of such an expression. In 
other words, we were not able to compute the next to 
leading order terms of this expression. 



< 1 
M 



£2 £3 
-M 



£1 £2 £3 




J gdg JM{tl0)JM{t20)Mh0) 

f ede M£ie)Jo(£29)J (£ 3 o) 



p 



1/2 



(cos 712 



P_' 1/2 (C0S71 2 ) 



cos Af7 12 (96) 



where 712 is the angle between Li and L2 if (L±, L2, L 3 ) 
forms a triangle. Such a limit is valid again for finite 
values of M only. 

It is then possible to transform the discrete sum on £ 3 
in Eq. (|71[) into a continuous integral on £3. From the 
expression of the Ponzano-Rcgge limit 



D. Examples 

We explore the consequences of our flat-sky expression 
of the bispectrum on a series of example in order to show 
that it is a robust expression in most practical cases. 



1. Flat case 

We first consider the so-called flat case for which the 
primordial bispectrum takes the form 



B aat (k 1 ,k 2 ,k 3 



(102) 



ll «2 «3 





1 



2-KA{L 1 ,L 2 ,L i ) 



(97) 



where A(Li, L2, L 3 ) is the area of the triangle formed by 
L\. L2 and L 3 , we can easily obtain that [26| 



2L, 



£1 
M 



£2 £3 
-M 



£1 £2 £3 




2tt 



>s(M/3 {L} ) b tlt2 t 3 



where /3{l} is the angle formed by Li and L 2 , i.e. 
Ll = L\ + Ll + 2L 1 L 2 cos/3 {l} . 



(98) 



(99) 



It has to be emphasized that this continuous limit can 
only be taken when not only £\ and £2 are large but also 
when their difference \£\ — £2] is large, so that the sum is 
not dominated by discrete values when /3{l} — > 7r. Given 
this limitation one gets by identification the expression 



of the bp 



in the flat-sky limit 



(100) 



In the Sachs- Wolfe limit, the bispectrum of the CMB 
temperature can be computed exactly and is given plj 
by 



Lflat 



L\ + L 2 + L3 



{L 1 +L 2 + L 3 ) 2 -9/4 



L1L2L3 

whereas our flat-sky limit (|100[) gives 

.fiat 2 



L\Li2L 3 {L\ + L2 + £3) 



(103) 



(104) 



These two expressions differ only by a term of the order 
of 1/(L 1 +L 2 + L 3 ) 2 . 

The computation of £^ 3 a/ is more complicated since 
there is no simple expression available. In the flat-sky 
approximation we get 



t-flat 



Mir — cmE 



■U1I3 _ 



d M K 



( Uxls \ 



(105) 

where E and K are Elliptic functions of the first and 
second kinds respectively, cm and dy[ are homogeneous 
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100 



FIG. 6: Left: £f^ M for £\ = 10, M = 0, 1, 2 (top to bottom) from exact computations (dots) and fiat-sky approximation (solid 
line) in the Sachs- Wolfe plateau limit. Right: ^°5m f° r ^i = 10, M = 0, 1, 2, 3 (from thick to thin dots) computed exactly from 
Eq. Gill- 



functions of £s/£i that are such that the numerator in with 



Eq. ([T05]) scales like (4/4 



\-l-M 



when t 3 j£\ is large, 



Co 

Cl 
C2 

<*2 



1 = -do 
1 = di 

(2l 3 /li + l)(4/li + 2) 
4/4 

(2(£ 3 /£ 1 ) 2 - 3£ 3 /h + 2) 



(106) 
(107) 

(108) 
(109) 



The result formally diverges for £\ — > £3, due to the fact 
that some configurations are IR divergent when k% — > 0. 
It implies that for the flat case in the Sachs- Wolfe plateau 
limit, the set of M's contributing in the sum (|8ip remains 
finite and this validates the flat-sky approximation for the 
bispectrum given by Eq. (jlOOp . The result is depicted on 
Fig.! 



2. Local case 



In the local case, the bispectrum is defined [2l| by 



D 



local 



(fci,fc 2 ,fc 3 ) 



1 



1 



1 



1,31,3 



ft 3 «, 1 



3 ■ 



(110) 



The bispectra can be computed by splitting B in 3 terms 
none of which has pathological IR divergences (in other 
words, bi x i 2 i 3 can be symmetrized only at the end) and 
we shall not encounter the divergence problems we had 
in the flat case. In its thin-shell approximation, a similar 
method consisting in splitting the computation in three 
terms was also of particular interest for the numerical 
computation of the bispectrum generated by non-linear 
effects around the LSS that we performed in Ref. Q. 

Furthermore, in the Sachs- Wolfe regime, the only non- 
zero contribution for each term is for M = 0. This leads 
to the flat-sky limit of the bispectrum 



[local ( ^local i ^local 

°e 1 e 2 e 3 — ls^.,0 + 



/^local \ 
' S£ 3 £ 2 0j 



(111) 



^local 



(112) 



Had we chosen the flat constraint L = k±r instead of 
L = fcj_r, we would have obtained a similar expression 
with L — > L. It would thus have matched the full-sky ex- 
pression, similarly to what was obtained in Section III CI 
Again, we conclude that for the local case, in the Sachs- 
Wolfe plateau limit, only M — is contributing and this 
validates the assumption that the number of M's con- 
tributing in the sum (|81|) is finite. Beyond the Sachs- 
Wolfe regime we can estimate numerically the different 
terms in the sum on M in order to validate the assump- 
tion that only a finite number of M is required to es- 
timate the bispectrum out of the t;e 1 i 2 M- Similarly to 
the Sachs- Wolfe plateau, we split in three terms the pri- 
mordial bispectrum (|110[) . and for each of these terms 
we plot ^jaif on Fig. |6] for different values of M. It 
is clear that the contribution of each M is exponentially 
supressed when M increases. 

To finish, we compare in Fig. [7] the flat-sky limit to 
the exact full-sky calculation for the bispectrum of local 
type and for equilateral configurations (i.e. such that 
£l = £ 2 = £3 = £). The agreement is excellent. 



VII. CONCLUSION 

This article provides a systematic construction of the 
flat-sky approximation of the angular power spectrum 
both for the temperature and polarization, in particular 
it shows that the expansion can be performed to any 
order. Additionally, we showed that this construction is 
not unique and that there exists a two-parameter family 
of flat-sky expansions, depending on the arbitrary choice 
of the flat-sky directions with respect to the azimuthal 
angle. 
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FIG. 7: Comparison of the flat-sky approximation (dashed 
line) to the exact full sky computation (solid line) for the 
reduced bispectrum but in the case of local type primordial 
non-Gaussianity. 



As long as the sources are scalar, the first correction 
term scales as l/£ 2 for a proper choice of flat-sky con- 
straints (k±r = £+1/2 or k±r = \J l{l + 1)). In more re- 
alistic cases involving direction-dependent sources, such 
as the Doppler term or the anisotropic stress, the cor- 
rections terms were shown to scale as l/£ whatever the 
flat-sky constraint and are given in Eq. (|53|) . 

Two particular extensions of the flat-sky expansion are 
particularly useful: the thin-shell and the Limber ap- 
proximations, depending on the spatial extension of the 
sources. We checked that the corrective terms obtained 
in the Limber approximation are consistent with existing 
literature Q , and we recovered the expression derived at 
leading order in the thin-shell approximation [l4| where 
all flat-sky expressions are identical. We discussed here 
the validity of those approximations. 

For practical purposes, the "best" formula to use de- 
pends on the context: 

1. for the CMB, as long as the thin-shell approxima- 
tion is a good approximation, that is on the LSS, 
one should use Eq. (f5U| since one cannot neglect 
the Doppler term. It includes corrections of order 
l/£, and its validity is limited by the errors intro- 
duced while taking the thin-shell limit which are 
at most of order AfLSs/^LSS < 0.01. They appear 
in practice to be less on small scales and on large 
scales they can be corrected with Eq. (|46|) . 

2. for large scale structures such as galaxy catalogs or 
weak-lensing, for which the Limber approximation 
is a good approximation, one should use Eq. flU]). 
It includes corrections of order l/£ 2 and its validity 
is limited to l/£ 3 [27]], as long as there are no larger 
errors introduced while taking the Limber approx- 
imation. Note also that for the late integrated ef- 
fects and the effect of reionization on the CMB, we 



should also use a Limber approximation. 

Generalization of this construction scheme to the bis- 
pectra was found to be more cumbersome. It led us to 
introduce an alternative description of the bispectra for 
which the flat-sky approximation is well controlled. This 
corresponds to the coefficients (U^t^M defined in Eq. (fTTj) 
whose relation with the usual reduced bispectrum form 
bg 1 g 2 e 3 can be found in Eq. (fTTj) and inverted in Eq. (|72p . 
For this quantity we were able to propose a well con- 
trolled flat sky approximation. It actually leads to a 
specific form of flat-sky approximation for b^ 2 ^ 3 , as de- 
scribed in Eq. in the sense given by Eqs. (jlOOlllOll) . 
the next-to-leading order terms of which remain however 
obscure (we encounter here exactly the same difficulty 
as when one tries to use the correspondance of Eq. (p}). 
In this case also one can further simplify the numerical 
integrations by using the thin-shell or Limber approxi- 
mations. 

The validity of the bispectrum flat-sky leading order 
expansion was tested numerically in simple cases of sep- 
arable primordial bispectrum such as the local and flat 
primordial bispectra. It was found to be very accurate 
in those cases, below the 1% level. Such expressions are 
obviously of great interest for non-separable shapes of 
primordial non-Gaussianity since no fast full-sky method 
is known yet. 
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Appendix A: Useful properties of the Bessel 
functions 

The Bessel functions satisfy the recursion relations for 
n > 1 

Ti 

J n (x) = J n -\{x) Jn(x), (Al) 

X 

2 ?i 

J„+i(x) = — J n (x) - J n -x{x) , (A2) 

with Jq = —J\. 

We define the integrals 

/>oo 

I n:P {a,b)= / x ( - 2p+n+1 '>J n (ax)Jo{bx)dx. (A3) 
Jo 

By taking successive derivatives with respect to a we ob- 
tain 



h tP (a,b) = -d a I , p (a,b) 



(A4) 



I 0tP+1 (a,b) = -d a [ah tP (a,b)]. (A5) 



19 



Using the orthonormality relation 



i S((i 6) 
xJo(ax) Jo(bx)dx = , (A6) 



we obtain the expression of the main integrals of interest 
for this paper 



h,\{a,b) 



6(b - a) 
b 

S'(b-a) 



5'{b-a) S"(b-a) 

ab b 
S'(b-a) 8"{b-a) 8"'{b-a) 



ba 2 



ab 



(A7) 
(A8) 
(A9) 
,(A10) 



Appendix B: Relating the corrections in the 
different flat-sky expansions 



In the expansion ([22]), had we chosen the constraint 
kj_R = L, then the expression of the corrective terms in 
the square brackets would have been modified by the re- 
placement of L by L and an extra operator — l/(8L 2 )Z?j_ 
would have been added in the brackets since for any func- 
tion Q, Q{L/r) ~ Q(L/r) - l/{8L 2 )kQ'(k). This means 
that the numerical factor in front of the operator D± 
would be 1/6 instead of 1/24. For a scale- invariant power 
spectrum and a sharply peaked constant transfer func- 
tion, i.e. w(k, r) oc S(r — tlss), the correction terms are 
thus proportional to 



(AD ± -D S ± ) 



1 



15k 4 , 



(fcj - 6fc r 2 ) 
fc 9 



and the integral on k r of this function is zero. This is a 
good check of our expansion. Indeed, in Section III CI we 
have seen that in this particular case, the leading order of 
the flat-sky expansion is equal to the exact result so that 
it was expected that the corrections vanish. Actually, this 
property should not depend on v since the w-dependence 
in the lowest order term of the expansion coming from 
the fact that R is a function of v, disappears when the 
sources are located only at a given distance tlss- For 
all choices of v and for a scale invariant power spectrum 
and a peaked transfer function, the lowest order of the 
flat-sky expansion is equal to the exact result. We can 
check indeed that in that case the corrective terms in 
the expression (|19|) where v = once expressed with the 
constraint k±R = L, vanish as well. The ^-dependent 
corrections need to be consistent, and though they are 
formally different as v is varied, they should lead to the 
same result, at least in that case where the transfer func- 
tion is sharply peaked. We now show that the corrective 
terms for different choices of v can be related by integra- 
tion by parts in that case. 

We assume that w(k,r) cx 5(r — ?*lss) so that we can 
also assume that all the k dependence is located in P(k) 



so that we can drop the factors if(fc,rLSs) f° r simplicity 
and set r^ss = 1- Expanding the correlation ([7]) in Z we 
obtain 



Cu - £.1/2 + 



dk, 

W) 

(1 - 2u)Z 3 



2 fcidfciP(fc) 
[-Zk 2 J a {k ± Z) + k ± Ji(k ± Z 



(Bl) 



and thus we can relate the corresponding multipoles 



C°, = d v - 1/2) + ACf 



with 



AC f = 



(l-2v) 



dk r 



dk ± P(k) 



(2tt 

h$_Ii,i(kx,L)-kxh£l , 2 (kx,L) 



(B2) 
(B3) 



(1 - 2v) 



dk r 



(2tt 

x lk 2 ± I hl (k ± ,L)-k 2 



dfc ± P(fc) 
d 



dk ± 



Integrating by parts we obtain 



ACf 



(1 - 2v) 



dk r 



dk. 



(B4) 



xkl { P(fc)ii,i(* ±J i) + ^P'(fc)/i,i(fc ± ,L) 



Another integration by parts in k r is sufficient to obtain 
that AC i = 0. This proves that even if the expression 
of the corrections is formally different for a different v, 
the corrective term remains the same in the thin-shell 
approximation. 



Appendix C: Asymptotic forms for the spherical 
harmonics 



The spherical harmonics can be expressed in terms of 
the Legendrc polynomials as 



rrf 


~{l-m)\~ 


1/2 




_(e + m)\_ 





(Cl) 

where u = cos 8. The construction of the spherical har- 
monics for non-zero values of m are obtained from the 
recursion relation 

d m 

pr(u) = (-ir(i - uy^—PM . (C2) 

Each of these operators can then be evaluated in the flat- 
sky limit, with u ps 1 — 9 2 /2, 



-V 

du J 



\9deJ 

m(m — 1) 

6 

6 



(C3) 



d 

dd0 
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Taking advantage that for finite values of m 



spherical harmonics, 



{l-m)\ 



it 



m) 



1/2 



1 



(£ + 1/2)" 
m(2m + l)(2m 
24£- 



(C4) 



we can use the relation (|T5|) to get the expression of Yj, m , 
and its next to leading order corrections to an arbitrary 
order, in the flat-sky approximation. The final result can 
be written in terms of Bessel functions of the first kind 
thanks to the relation, 

(j \ m 
m) Md) - (C5) 

We then recover the known large £ expression of the 



2L 



1/2 



xrvM ->• ( ^ ) (-iy n J m [L9w m ^ (ce) 



where the first correction term is in l/£ 2 or in 9 2 . It 
comes from the sub- leading terms in Eq. (TT5)) . the sub- 
leading operators in (|C4|) and the subleading terms in 
(|C4[) . All these corrective terms make the flat-sky ex- 
pansions controllable when m is finite but clearly not 
when it is of the order of I. The corrective terms are 
expected to involve Bessel functions J m i where ml differs 
to m by at most 2. Subleading terms in the final results 
of flat-sky expressions can then be obtained in principle 
using relations (or similar to those) given in appendix [A"l 
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